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i t Active pre-main-sequence binaries with separations of around ten stellar radii 

present a wealth of phenomena unobserved in common systems. The study of these 
objects is extended from Classical T Tauri stars to the Herbig Ae star HD 104237. The 
primary has a mass 2.2 ± 0.2 M and secondary 1.4 ± 0.3 M Q . Spectro-interferometry 
with the VLTI/ AMBER in the if-band continuum and the Br7 line is presented. 
It is found that the if-band continuum squared visibilities are compatible with a 
, circumbinary disc with a radius of ~ 0.5 AU. However, a significant fraction (~ 50 per 

cent) of the flux is unresolved and not fully accounted by the stellar photospheres. 
The stars probably don't hold circumstellar discs, in addition to the circumbinary 
disk, due to the combined effects of inner magnctosphcric truncation and outer tidal 
truncation. This unresolved flux likely arises in compact structures inside the tidally 
disrupted circumbinary disc. 

Most (> 90 per cent) of the Br7 line emission is unresolved. The line-to-continuum 
spectro-astrometry shifts in time, along the direction of the Lya jet known to be driven 
by the system. The shift is anti-correlated with the Br7 equivalent width. It is shown 
that the unresolved Br7 emission cannot originate in the jet but instead is compatible 
with stellar emission from the orbiting binary components. The increase in the absolute 
value of the equivalent width of the line takes place at periastron passage; it could arise 
in an accretion burst, a flare or in the increase in effective size of the emission region 
\ by the interaction of the magnetospheres. The binary longitude of the ascending node 

is found to be ft — (235 ± 3)° and the orbit retrograde. 

The origin of the jet is revisited. The tidal disruption of the circumstellar disks 
creates difficulties to ejection models that rely on stellar magnetosphere and disc cou- 
pling. A scenario of a stellar wind collimated by a circumbinary disc wind is suggested. 
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1 INTRODUCTION 

Pre-main-sequence actively accreting binaries with separa- 
tions of tens of stellar radii present a wealth of phenomena 
such as tidally disrupted discs, colliding magnetospheres or 
periodic accretion. In parallel, these objects present overall 
properties (such as line profiles) that are very similar to their 
single counterparts. This is puzzling given that the close en- 
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vironment is severely disrupted by the orbital movement of 
the system. The evidence for tidally disr upted circumbinary 
actively accreting discs is growing (e.g . iJensen fc Mathieul 
1 19971 ; lAndrews etai1l201ll ; Harris et al.ll2012l ). but for only 
very few objects there is enough information allowing the 
study of the effects on the several components of the pre- 
main-sequence star paradigm such as the magnetosphere, 
the disc, the winds and jets. In this article, this sample is 
extended by presenting a study of HD 104237. 

The HD 104237 spectroscopic binary has parame- 
ters P = 20 d, a = 0.22 AU, mi = 2.2 M 



Bohm et al. 2004, Ap pend ix |A")> . a luminosity L ~ 35 Lq 
1l 19981) and is loca ted at a distance of 



van den Anckcr et al. 



d = (116 ± 7) pc ( Perrvman et al.lll997l) . It presents strong 
varia ble emission in Ha (jBaines et a l. 2006; Bohm et al.l 
l2006h . Using the Br7 line. lGarcia-Lopez et all l|2006h derived 
an accretion rate of M acc ~4x 10 -8 Mq yr _1 . The spectral 
energy distribution (SED) was fitted with a circumbinar y 
ring with an inner rim at 0.45 AU bv lTatulli et al. | i|2007l ). 
This source also drives a l arge-scale jet detec ted in Lya at 
position angle jc t = 332° (| Grady et al.ll2004l ). 

Most of the Br7 emission from young stars is thought 
to ar ise in the stella r magnetosphere (e .g. iKurosawa et al.l 
|2011| ). The report bv lTatulli etaD (|2007l ) of HD 104237 Br7 
originating at sc ales as large as th e inner rim was surpris- 
ing. In contrast, iKraus et alj (|2008h found that the Br7 for 
HD 104237 originates in a region smaller tha n the contin- 
uum. In the survey of Br7 in young stars of lEisner et al.l 
l|2010l ) most objects have line emission more compact than 
the continuum emitted by the disc rim. However, some ob- 
jects do present a contribution from ex tended Br7, for ex- 
ample, more massi ve Herbig Be stars (jMalbet et al.ll2007l ; 
IWeigelt et aill201ll ). The previous HD 104237 observations 
had a poor (u, v) coverage; the unavailability of fringe track- 
ing resulted in low signal-to-noise ratio (SNR). This partial 
information could only approximately constrain the origin 
of the Br7 emission. Is it originating in the magnetosphere, 
the gaseous disc, the inner rim or jet? Is HD 104237 Br7 
emission somehow related to the increase d activity in tight 
systems as found in UZ Tau E and DQ T au l|Basri et al.ll 19971 ; 
I Jensen et al] 1200/1 : iKospal et aTll201lh ? These are some of 
the questions addressed in the present article. 

In Section [5] the observations and data reduction are 
presented. Then we proceed with the continuum emission 
results in Section [3] followed by the Br7 results in Section [4] 
The discussion starts with the Br7 origin in Section [S] and 
then concludes on the overall system in Section [6] Details on 
the stellar parameters used in this work and the data analy- 
sis method used to extract line-to-continuum visibilities are 
presente d in the App e ndix. Throughout the paper the no- 
tation of iGlindemannl |201ll) is used for the interferometric 
quantities. 



2 OBSERVATIONS AND DATA REDUCTION 

The observations were conducted at the VLTI, using the 
near- infrared instrument AMBER l|Petrov et al.l 120071 '). in 
the spectral medium-resolution mode (MR; R ~ 1500) and 
the spectral high-resolution mode (HR; R ~ 10000), cover- 
ing the Br7 line. The MR observations were performed with 
the auxiliary telescopes (ATs), the HR observations with the 



unit telescopes (UTs). The field of view is about 250 mas 
for the ATs and 57 mas for the UTs. The FINITO fringe 
tracker (|Le Bouquin et al"1l2008h was used throughout the 
observations, although it was operated at the limit of its sen- 
sitivity with the ATs. To calibrate the continuum, we also 
used low-resolution data (LR; R ~ 35) obtained without 
fringe tracking and with the ATs. The baseline lengths var- 
ied between 14.5 to 125.3 m, corresponding to a maximum 
angular resolution \/2B of 1.8 mas. Each measurement was 
interleaved by observations of a calibrator star ( HD 118934, 
diameter 0.877 ± 0.005 mas, iMerand et ail 120051 ). Each MR 
and HR data-set typically consists of 5 to 10 individual files, 
each with roughly 50 individual exposures (or frames). The 
LR data-sets consist of a sequence of 5 files of 1000 frames 
each. 

The data reduction was performed using the AMDLIB 
package (v3.0b), fol lowin g standard p r ocedu res described in 
iTatulli et ail (|2007l ) and IChelli et all (|2009l ). Raw spectral 
visibilities, differential phases, and closure phases were ex- 
tracted for all frames of each observing file. For each data 
set and baseline, the distributions of the fringes SNR ver- 
sus optical path delay, the visibility versus SNR and his- 
tograms of the SNR were inspected. For the MR and HR 
data, frames with the same instrumental setup and close in 
time were merged to enhance the SNR. For all data-sets, a 
selection of 20 per cent of the highest quality (SNR) frames 
was made. Other selections were done without appreciable 
differences in final results. The transfer function was ob- 
tained by averaging the calibrator measurements, after cor- 
recting for its intrinsic diameter. In practice, we found that 
the differential phase root-mean-square (rms), rather than 
the visibility rms, is a good indicator of the data quality. 
For a well-resolved object, the visibility rms will be small 
but not indicative of the data quality. We therefore applied 
a second criterion for data selection: only data with a rms 
in the differential phase smaller than 27r/50 rad was kept. 
After data reduction, we used a custom so ftware to calibrate 
the wavelength grid (jMerand et al.ll2010r ). The accuracy of 
the wavelength/velocity is < 50 km s _1 . A summary of the 
retained data is presented in Table [1] The UV coverage of 
the final data set is presented in FigfT] 

The absolute value of the visibilities obtained with the 
UT baselines could not be determined due to random vibra- 
tions of the telescopes and/or different performance of the 
fringe tracker during the target and calibrator observations. 
We also found some discrepancies between the absolute val- 
ues in the continuum visibilities retrieved from the MR/ AT 
observations and from the LR/AT data. This issue affects 
all spectral channels in the same way and does not change 
any conclusion relative to the spectral analysis. The result of 
the MR/HR data reduction is therefore squared visibilities 
normalized to the continuum (Vr/c = V^/Vq) and differen- 
tial phases relative to the continuum (A0 T /c — 0t — 0c)- 
The final data-set is presented in Fig.[2]for the MR/ AT data 
and Fig. [3] for the HR/UT data. 

In Fig. [4] the Br7 spectra are presented. The data for 
each night was continuum normalized and averaged. The 
first three nights spectra (MR1, MR2 and MR3) were ob- 
tained with medium spectral resolution, the fourth night 
spectrum (HR1) with high spectral resolution. The line in- 
tensity increases with time but the profile shape remains 
very similar. 
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Table 1. Log of the observations. 



Spectral 


Date 


Telescope 


B (ra) 


PA (°) 


Spectral 


Date 


Telescope 


B (m) 


PA f°1 


Mode 




Stations 






Mode 




Stations 






HR 


2010 Mar. 04 (HR1) 


U2-U3 


35.4 


70.0 


LR 


2008 June 03 (LR2) 


H0-G0 


29.6 


-69.8 






U3-U4 


60.5 


115.6 






G0-E0 


14.9 


—69.8 


MR 


2010 Feb. 07 (MR1) 


G1-D0 


57.1 


-72.4 






H0-E0 


44.4 


-69.8 




2010 Feb. 10 (MR2) 


K0-G1 


64.7 


-155.4 


LR 


2009 Apr. 04 (LR3) 


G1-D0 


61.2 


-58.6 






G1-A0 


82.4 


-91.6 






D0-H0 


62.6 


68.5 






K0-A0 


125.3 


-119.2 






G1-H0 


55.1 


6.2 




2010 Feb. 12 (MR3) 


H0-G0 


30.6 


-93.5 


LR 


2009 May 21 (LR4) 


G1-D0 


62.6 


-27.6 






G0-E0 


15.3 


-93.5 






D0-H0 


58.3 


54.0 


LR 


2008 May 26 (LR1) 


G1-D0 


59.9 


-37.1 






G1-H0 


40.0 


24.5 






D0-H0 


56.2 


90.0 


LR 


2009 May 22 (LR5) 


HO-DO 


57.1 


-56.4 






G1-H0 


49.3 


24.2 






DO-AO 
HO-AO 


28.5 
85.6 


-56.4 
-56.4 



B and PA are the average baseline length and position angle. In the text, the nights are referred to as night HR1 (2010 Mar. 04), night 
MR1 (2010 Feb. 07), night MR2 (2010 Feb. 10), night MR3 (2010 Feb. 12), night LR1 (2008 May 26), night LR2 (2008 June 03), night 
LR3 (2009 Apr. 04), night LR4 (2009 May 21) and night LR5 (2009 May 22). 
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Figure 2. MR velocity-corrected data (black) and best model fit (red, cf. Section|4]and Appcndix[B}. The rose background is the region 
used in the fit, the blue legend contains the night code, baseline and hh:mm of the data. Top: continuum-normalized squared visibility 
(Vr^y c ); bottom: differential phase (A</> t /q). Appendix IB3I discusses the fit quality and in particular the data set MR3 H0-G0 07:51. 
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Figure 3. HR velocity-corrected data (black) and best model fit (red, cf. Section [4] and Appendix [Bj . The rose background is the region 
used in the fit, the blue legend contains the night code, baseline and hh:mm of the data. Top: continuum-normalized squared visibility 
(Kfy c ); bottom: differential phase (A0x/c)- 
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Figure 1. UV coverage of the observations, LR data (red, •), MR 
data (black, •) and HR data (blue, •). The dashed line position 
angle depicts the large-scale jet direction (0j c t = 332°). 
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Figure 4. Spectroscopic observations for each night. Medium- 
resolution: night MR1 (red), night MR2 (green), night MR3 (blue) 
and high-resolution: night HR1 (magenta). 



3 THE CONTINUUM EMISSION OF HD 104237 

The 7f-band continuum data consists of LR squared visibil- 
ities (V 2 ) and LR and MR closure phases. The MR and HR 
V 2 continuum values are unavailable due to transfer function 
variability. Only one set of MR data allows the computation 
of the closure phase. The data is presented in Fig. [5] In the 
following the best fit geometric model to the data is pre- 
sented, then the scatter in the data is explained with the 
binarity of the object. 

Geometric model In the astrophysical context of 
HD 104237, the V 2 continuum data was fitted with an ex- 
tended ring concentric with a point source (Model CI). The 



Table 2. Continuum V 2 models and derived parameters. 



Model Value (mas) / (per cent) \ 2 

Model CI r = 3.5 ± 0.8, a = 1.8 ± 0.6 45 ± 7 5.4 
Model C2 r = 4.0, a = 1.8 60 n.a. 



Models CI and C2: ring of radius r with radial Gaussian <r 
cross-section and point source with fractional flux f;xt ' ls the 
reduced \ 2 ■ 

extended ring has a radius r and radial Gaussian ring cross- 
section. It mimics an unresolved spectroscopic binary sur- 
rounded by an extended circumbinary sublimation rim. The 
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Figure 5. if -band continuum data. Top: V 2 versus spatial fre- 
quency radius R in the (u, v) plane. The Model CI line is the best 
fit to the data. The rose range and Model C2 depict the effect of 
the binary (cf. Section[3]for details). Bottom: closure phase (CP) 
versus maximum triangle spatial frequency radius (i^max 

). The 

data colour codes are: night LR1 (red, ■), night LR2 (green, +), 
night LR3 (blue, A), night LR4 (cyan, •), nig ht LR5 (mag enta, 
x), night MR2 (black, • in closure phase onlv). lKraus et alj data 
is plotted in orange. 

best-fitting parameters and reduced x 2 ar e presented in Ta- 
ble [2] and the data and fit in Fig. [5] The fit presents a sig- 
nificant xt an d should be taken as a rough estimate of the 
circumbinary disc size and unresolved flux fraction. The esti- 
mated disk size from model fitting is model dependent. Using 
other models we found it to vary from 3.0 to 4.4 mas. A re- 
sult common to all models is that the fractional flux found in 
the point component is larger th an the 30 per c ent ex pected 
from the SED fits to the disc by ITatulli et all (|2007h . 

Effect of the spectroscopic binary Model CI predicts 
zero closure phase. The binary asymmetric emission gener- 
ates a non-zero closure phase. To explain the closure phase 
the unresolved flux of Model CI was allocated to both com- 
ponents of the binary (described in Table IA1|I . Fitting the 
closure phase with such a model had 3 < xt 5, which sup- 
ports the spectroscopic binary as causing most of the ob- 
served closure phase values and scatter. 
With regards to the V 2 , the effect of including a binary is 
to decrease the visibility. The projected time averaged bi- 
nary separation in the plane of the sky is (2.2 ± 0.7) mas. 
The binary is marginally resolved at the longer baselines. In 



Fig. [5] (Top) the effect of including a binary in the modeling 
is illustrated. Model C2 has an increased radius r = 4.0 mas 
and increased unresolved flux fraction / = 60 per cent. Ap- 
parently this model does not fit the data. However, if the 
unresolved flux of Model C2 is allocated to the binary then 
the resulting range of V 2 values (rose range in the Figure) 
covers most of the data points. The binary is a satisfactory 
explanation to the observed scatter in the V 2 data. 

Summary of the fits The continuum interferometric data 
of HD 104237 is reasonably well explained by a spectro- 
scopic binary surrounded by a circumbinary disc of radius 
~ 4 mas. This v a lue is intermediate to the 3.9 mas derived 
by iTatulli et alj (|2007l ) and 5.0 mas found by iKraus et al.l 
(2008). The flux allocated to the circumbinary disc is smaller 
(/ ~ 50 per cent) than the if - band excess ( / ~ 70 per cent) 
derived from the SED fits by ITatulli etalf (|2007l ). Signifi- 
cant excess exists in unresolved emission in all models. The 
continuum emission will be discussed in Section \6. II 



4 THE BR 7 COMPLEX VISIBILITIES 

4.1 Line-to-continuum observables 

As the SNR of the data is relatively low and the line 
marginally spectrally resolved in MR, a special procedure 
was developed to extract the line-to-continuum visibilities 
and phases. This procedure called least-squares deconvolu- 
tion is presented in detail in Appendix [B] It provides a fac- 
tor of two increase in the precision of the extracted line- 
to-continuum squared visibility and differential phase. In 
Fig.[S]the extracted quantities are presented. The overall re- 
sults are that the Br 7 is more compact than the continuum 
Vl/c > 1 an d that the line centroid presents sometimes a 
displacement with respect to the continuum A0 L /ci 7^ 0. The 
line-to-continuum closure phase for the MR2 night, where a 
complete triangle was available, is found to be zero within 
errors CP = (0 ± 7)°, i.e. in the line-to-continuum spatial 
asymmetry ("skewness") is the same. 

4.2 Spectro-astrometry 

The shift property of the Fourier transform relates the 
measured wavelength-differential phase to the spectro- 
astrometric centre of the source brightness distribution. In 
particular, the line-to-continuum phase A0 L /c: is related to 
the line-to-continuum spectro-astrometry £*l/c v i a 

= -27TQ! L/C ■ R, 

with R the position vector in the (u, v) plane. HD 104237 is a 
spectroscopic binary with a period of 20 d. Therefore, fits for 
each night are presented in Fig. [7] (top). Except for the lowest 
angular resolution data (smallest baselines), the measured 
shifts are displaced from the origin by more than 3er. Because 
these are line-to-continuum positions the observed shift can 
be caused either by a continuum shift, or by a line shift or by 
both. The reader is reminded that in Section [3] it was found 
that the continuum had a non-zero closure phase, therefore 
its astrometric centre is not at the origin. 

Fig. (top) shows that the spectro-astrometric position 
changes from night to night along the Lya jet direction. 
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Figure 6. Extracted quantities from the least-squares deconvo- 
lution for each baseline versus frequency radius R in the (it, v) 
plane. Line-to-continuum visibility ratio Vl/c (top) and line-to- 
continuum differential phase A<^>l/c (bottom). Note that Vl/c i s 
a ratio and can be larger than 1. MR data of night MR1 (red, 
■ ), night MR2 (green, +) and night MR3 (blue, A). HR data of 
night HR1 (magenta, •). 



Figure 7. Br7 spectro-astrometry. The continuum is assumed 
at the origin. Top: average position (relative declination A<5l/c 
and right ascension ARA L /c) f° r an nights' data (black, •fc); 
for the individual nights the coding is the same as in Fig. [6] The 
ellipses present the lc error in the spectro-astrometry. The dashed 
straight line shows the jet position angle. Bottom: shift along jet 
axis versus Br7 equivalent width. 



Comparison with the integrated spectra (Fig. [7] bottom), 
which are also variable with time (cf. Fig. [4} shows that the 
shift magnitude correlates with the Br7 equivalent width. 
Larger Br7-to-contimium flux translates in smaller shifts. 

4.3 Complex visibility 

The line-to-continuum wavelength differential phase and vis- 
ibility can be fitted simultaneously. The model for the com- 
plex visibility is 

Vl /c = exp(-2™* L /c • *) y o( ^ cl y 

It consists on a displacement cvl/c an d a centre-symmetric 
visibility Vl whose functional form is described in Table [3] 
Model CI is used to calibrate the data, as it is centre- 
symmetric it has a zero closure phase. This model is not 
necessarily optimal because, as shown in Section [3] the con- 
tinuum has a non-zero closure phase. However, a calibration 
of the line closure-phases in the same way as was done for 
the visibility would add too many parameters and make the 
result even more model dependent. For a given night only a 
small number of (u, v) points is available and there is poor 
baseline distribution. Therefore, a night-by-night fit of the 



data was not attempted. With this caveat in mind, the fits 
are presented in Table |3] and Fig. [8] The main result is that 
the Br7 emission is composed of an unresolved component 
accounting for at least of ~ 90 per cent of the flux. Given the 
uncertainties in the continuum visibility calibration, a model 
where all the Br7 flux is unresolved (Model M3) is also 
roughly compatible with the data. The size of the extended 
remaining Br7 component is compatible with the contin- 
uum fit presented in Section [3] but this result is strongly 
dependent on the actual geometric model used. We there- 
fore abstain from over-interpreting this "resolved" compo- 
nent of Br7 obtained in the models of Table [3] The models 
also present an astrometric signature at the position of the 
average spectro-astrometry. 

4.4 Summary of the modeling 

The results of the modeling can be summarized as follows: 

- The Br7-to-continuum spectro-astrometry position 
moves along the jet direction. 

- The Br7-to-continuum spectro-astrometry position 
changes with time. 

- The Br7-to-continuum spectro-astrometry movement 
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Table 3. Line-to-continuum complex visibility models and derived parameters. 



model 


rel. 8 (mas) 


rel. a (mas) 


(mas) 


/ (per cent) 




Model Ml 


-0.17 ±0.04 


0.38 ±0.08 


r = 5.4 ± 1.0 


92 ±1 


6.7 


Model M2 


-0.17 ±0.04 


0.38 ±0.08 


a = 3.5 ±1.6 


90 ±3 


6.9 


Model M3 


-0.17 ±0.06 


0.38 ±0.11 




100 


14 



Model Ml: displaced <5-ring of radius r and point; Model M2: displaced Gaussian and point source; Model M3: point source. / is the 
fraction of total flux in the point source. The continuum model used is Model CI from Table [2] 
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Figure 8. Line-to-continuum complex visibility models. The lines 
are the best-fitting models to the data. The fit information is pre- 
sented in Table [3] For A</>l/c the lines overlap. The data coding 
is the same as in Fig. \E\ 



is anti-correlated with the Br7 equivalent width absolute 
value. 

- Most (> 90 per cent) of the Br7 emission is unresolved. 

- The remaining extended Br7 flux is compatible in size 
with the continuum but this result is model dependent. 



5 THE ORIGIN OF THE BR 7 EMISSION 

In this section, we argue that the compact Br 7 emission 
observed along the direction of the large-scale jet cannot 
originate in the jet or a disc wind but instead arises in the 
close stellar environment of each binary component. Mixed 
models of binary plus jet were not attempted due to the high 
number of parameters with regard to the available data. 



5.1 Jet scenario 

The detection of Br 7 in HD 104237 along the direction 
of the Lya jet could possibly point to a common origin. 
Jets are known to emit hydrogen recombination lines, but 
the Br 7 em i ssion is normally c o mpact and unresolved (e.g. 
iDavis et all 1201 lT ). iBeck et all (|2010h find the Br 7 emis- 
sion to have a 2-10 per cent contribution from extended 
components at scales < 200 AU. At much smaller scales, 
Br 7 jets wer e detec ted in the 2008 Z CMa outburst by 
iBenistv et alj Jioiot). but most objects show unresolved 
emission (|Eisner et alll20ld : iKraus et al1l200St ). 



Model The Br 7 is modeled with a mixed contribution to 
the line from a displaced unresolved knot (Fknot) an d from 
variable stellar emission (F*) at a fixed position r*. The 
observed position (p) of the Br 7 is 



Ft , .Fknot 



with Ft 



F* + F k: 



and r kn 



dknotitjet directed 



along the jet position angle (9j ct = 332°). In Fig. [7] (top) 
the points MR1 and MR2 are located along the jet posi- 
tion angle, the point MR3 is compatible with a position 
at the origin. To explain this last point must be lo- 
cated opposite to the jet direction. The position angle of 
the star is therefore fixed at 6j ct — n, leaving free only d* 
(r* = — dirUjetj- The knot moves with speed u kn ot, its dis- 
tance is dknot = dknot,t + v knot{t — *<))• 

Visibility calibration The Br 7 visibilities and phases are 
recovered making use of the continuum Model CI. The line 
visibility (VL) then becomes 



Vl = VL/c Vc (Model CI) 

and the line phase ($>l) 

0l = 0l/c + 0c(ModelCl) = 



'L/C: 



(1) 



(2) 



because Model CI is centre symmetric 0c(ModelCl) = 0. 
This calibration has some limitations. Model CI overesti- 
mates the continuum squared visibility at small baselines 
(cf. Fig. [5} ■ The line visibility at small baselines is slightly 
overestimated. With regards the phase, the stellar compo- 
nents could be moving during the orbit. 

Fit The Br 7 visibility and phase are fitted with this model. 
The data and model are presented in Fig. [9] The best fit 
model parameters are presented in Table [4] The knot has a 
very high flux (Fk no t ~ 3.5 A), similar to values observed in 
a single emitting pre-main-sequence star. It is located very 
near the position of maximum astrometric displacement. It 
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Figure 9. Jet scenario. The night of observation is colour encoded as in Fig. [4] Line visibility (left) and line phase (centre) versus spatial 
frequency radius R in the (u, v) plane. The open squares are the best-fitting model. Right: drawing of the best-fitting model. The star 
(*) is the stellar position, the filled circle is the knot position, the crossed points the astrometric data. 



Table 4. Jet scenario best-fitting model parameters. 



Parameter 


Value 


F knot (A) 


3.46 ± 0.07 


dknot.to ( mas ) 


1.0 ±0.2 


«knot (mas d" 1 ) 


0.00 ±0.04 


(mas) 


1.0 ±0.6 




1.4 



has no proper-motion. The star is positioned at an almost 
symmetrical distance from the centre of the field-of-view. 

Physical inconsistency of the fit Although the model 
provides a reasonable fit to the data (xt = 1.4) it is phys- 
ically inconsistent as it requires a too high F^ not and too 
small Vknot- Firstly, Fk no t can be limited by an analysis 
of the line profile. As the system is seen almost pole-on, 
a significant flux in the knot would distort the line pro- 
file shape, with a blue-shifted enhanced contribution at 
the jet radial speed. Our analysis of the line profile limits 
i*knot < 1 A (12-25 per cent of the total Brr/ flux), a much 
smaller value than the best-fitting Fk no t • The Br7 line profile 
shape is simila r to those commonly observed in pre-main- 
sequence stars ( Folha fc EmersorJ l200ll ; iGarcia-Lopez et all 
l200fj ; iNaiita et al.1 119961 ). Secondly, in the best fit model 
the knot is at rest, at odds with the p roper-motion ex- 
pected for a jet knot. iGradv et al l (|2004h measure for the 
Lya knots HH 669- A and HH 669-B, a radial velocity of 
(—337 ± 80) km s _1 . The corresponding speed in the plane 
of the sky would be ui = 0.5(uj c t/360 km s _1 )mas d _1 . A 
third difficulty is related to the fixed stellar position. The 
central object is a spectroscopic binary, both the longitude 
of the ascending node and the time of periastron passage 
would have to be fine-tuned to keep the star in the required 
position opposite to the jet direction. 

We therefore rule out that a knot in the jet as the ex- 
planation for the compact Br7 emission. 



calibration to recover VL- The calibrated data is then fit- 
ted with orbital models in which the Br7 flux is allowed to 
vary. The resulting variable emission is then discussed and 
interpreted as variable accretion. 

5.2.1 Spectro-astrometry 

Assume that only the primary and the secondary contribute 
to the total Br7 emission and that the contribution is unre- 
solved. The primary fractional emission is defined as 

Br7i 

ql = — — 

Br7i + Br72 

and the primary fractional mass as 
mi 
mi + m2 

It can be showrQ that the Br7 spectro-astrometric centre 
(c) is given by 

C = (Ql — 0!m)fT/2, 

where T7/2 is the relative position of the primary with re- 
spect to the secondary. 

The previous equation for the centre is degenerate. The 
position angle of the orbit in the plane of the sky is deter- 
mined by the longitude of the ascending node (fi). Orbits 
with opposing Q. can provide the same astrometric centre 
provided the stars fractional fluxes are adjusted. Quanti- 
tatively, the values cil and Q yield the same astrometric 
centeio as the values a' L = 2a m — «l and Of = fi + 7T. This 
degeneracy is broken by the use of the visibility (VL)- Only 
in the very special case where 2a m = 1 is also the visibility 
symmetric. For HD 104237 2a m — 1.2T0.1 so both solutions 
are kept because they have similar xt (cf . Table . 



5.2.2 Visibility calibration 

As for the jet scenario we have recovered calibrated line vis- 
ibilities and phases making use of the continuum Model CI 



5.2 Binary scenario 

In this section we start by introducing the spectro- 
astrometric model of the data, then we present the Vl/c 



1 The spectro-astrometric centre is c = «L r i + (1 — Ch) r 2. 
The result is obtained by substituting the two-body solutions 
(7*2 = — a m r^ 2 and r\ = (1 — o m )ri/ 2 ) m the centre expression. 

2 As long as 2a m — 1 < «l < 2a m . 
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(cf. Eqs.[TJand[2]). This calibration does not address the fact 
that the continuum emission from the stars is moving during 
the orbital phase, creating a <j>c 7^ 0. It can be shown that 
this effect is 

cc = (011 — (011 + ct2)a m )r 1 / 2 = Qci'1/2, 

where a; = Fc,i/FT is the i th component continuum frac- 
tional flux and Ft = Fc,i +Fc,2 +-Pc,ciisc is the total contin- 
uum flux. The net effect of the continuum centroid drift due 
to the orbital motion of the stars is to bias ax with respect 
to the real one by roughly ag ■ HD 104237 has ai = 0.2 and 
q 2 = 0.1 (|Tatulli et alj|2007f ). thus a c ~ 0.02. A further 
fraction of ~ 0.2 of the total continuum is found unresolved 
in the continuum fits of Section O in worst case it would 
change the previous value to |qc] ~ 0.1. This bias is gener- 
ally smaller than our precision in ax and will only affect the 
fractional flux allocated to the binary components. 



5.2.3 Models 

The model has six parameters: two orbital (tp,i, Q) and four 
for the nightly variable flux ratio «l. The system orbit is 
computed using the parameters in Table IA1I The model is 
fitted to the Br7 visibilities and phases in two steps. In the 
first, the full x 2 space is searched for a global minimum. In 
the second, the minimum is fine-tuned with a Levenberg- 
Marquardt algorithm. Models with prograde and retrograde 
motion were fitted. The prograde models are a worse fit 
(xt = 4.0 and 4.8, for fl and Q.' , respectively) than the 
retrograde model (with xt = 3.2 and 3.4). The retrograde 
motion results are presented in Fig. [10] and Table [5] they are 
our baseline model for the interpretation of the Br7 emission 
from the system and are discussed in more detail in the next 
sections. 

The average time of periastron passage obtained in 
the above retrograde fit (tp,i) is independent of the time 
of periastron passage obtained in the radial velocity fit 
(tp, Appendix [A}. They are in very close agreement 
(tp — tp.i)/P = 2601.004. The spectro-astrometric data can 
be used to determine a more accurate longitude of the as- 
cending node by keeping tp,j fixed at the radial velocity value 
in the above fits, obtaining = (235 ± 3)°. 



5.2.4 Bry variability 

The total Br7 equivalent width (Br7T) derived from the 
AMBER spectra is variable (cf. Table [S] and Fig. 0). The 
interferometric data constrain the size (unresolved) and po- 
sition of the Br7 emission. The higher Br7 emission (nights 
MR3 and HR1) takes place near periastron passage (cf. Ta- 
ble [5|. The fit extracts the individual components emis- 
sion for each night which, as shown in Table O is variable. 
iDonehew fc Brittainl (|201ll ) refer to Br7 line variability in 
their sample of Herbig Ae/Be st ars of up to 4 A, when com - 
par ed with the observat ions of iGarcia- Lopez et al.l (2006) 
and lBrittain et all ({2007) . This equivalent width variability 
cannot be caused by the continuum. Roughly half of the 
continuum emission originates in the circumbinary disc in- 
ner rim whose Keplerian time-scales are ~ 70 d, much larger 
than the orbital period of the binary ~ 20 d. From night 
MR2 to night MR3 (48 h apart) a 3.5 A (90%) variation in 



the total Br7 equivalent w idth is measured. Furthermore, 
IDonehew fc Brittainl |201lf ) refer for their sample of Herbig 
Ae/Be stars a continuum variability of only < 10 per cent. 
The pole-on geometry of HD 104237 reduces shadowing ef- 
fects and continuum variability. The Br7 equivalent width 
variability is most probably intrinsic to the line and origi- 
nates in the binary stars very close environment. 



5.2.5 Accretion burst at periastron passage 

In classical TTauri stars the Br7 line could be formed 
in a magnetosphere. stellar wind or disc wind (e.g. 
Kuros awa et al.l 20 111 ; iKwan fc Fischerl 1201 il l. However, 



Kurosawa et al.l |201ll ') don't find signatures of the disc wind 



in the Br7 line, the bulk of the Br7 emission takes place 
on scales smaller than 0.1 AU, mainly from the magne- 
tosphere, and are unr esolved by the pre sent observations. 
This i s agreement with lKraus et all (|2008f ) and lEisner et al.l 
(2010) interferometric surveys of Br7 who also find com- 
pact emission compatible with magnetospheric accretion in 
classical TTauri stars and Herbig Ae stars. 

Although the paradigm of magnetospheric accre- 
tion is less solid for He rbig Ae stars (cf. discussion in 
Donehew fc Brittainll201ir ). a correlation between accretion 
rate and Br7 luminosit y was found from classical TTauri 
stars to Herbig Ae stars jMuzerplle et lii]|l99£j : ICalvet et al.l 
12004 IDonehew fc Brittairj|201ll ). It canbe argued that the 
Br7 variability is at l east an indirect trac e r of v ariable ac- 
cretion onto the stars. IGarcia- Lopez et alJ (|2006h measured 
for HD 104237 a Br7T value of ~ -5 A in 2004 and derived 
an accretion rate M acc ~4x 10~ 8 Mq yr" 1 . The exact frac- 
tion of the flux of each component could be slight biased by 
the continuum (cf. Section [5. 2. 2 [) . For the fit with f2 = 234°, 
the primary average Br7 flux more than doubles. For the fit 
with Q,' = 44° , it is the secondary Br7 emission that has an 
even stronger increase. There is therefore at least a twofold 
increase in the accretion rate from nights MR1 and MR2 to 
the nights MR3 and HR1 (as indicated by the spectra alone). 
The later nights are those in which the system is nearer peri- 
astron. Indeed, if we take only the MR data (which is taken 
over the course of 5 nights and is more tightly correlated), 
the system is approaching periastron and an accretion burst 
from night MR2 to night MR3 would take place. The night 
HR1 high accretion rate measurement takes place one period 
after MR3, again near periastron (Fig. llO|) . We cannot differ- 
entiate if there was a new burst of accretion at the second 
periastron passage or if the high accretion state remained 
during the full orbit. 



6 ACCRETION-EJECTION IN HD 104237 

This section presents a synoptic view. The tidally disrupted 
circumbinary disc is initially addressed, followed by the ac- 
cretion streamers and the interaction of the stellar magne- 
tospheres. Finally, the origin of the Lya jet is discussed. 



6.1 The tidally disrupted circumbinary disc 

The tidally disrupted inner AU is schematically represented 
in Fig. 1111 An inner gap is present and a circumbinary 
disc surrounds the system. The background image (from 
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Figure 10. Retrograde binary scenario. Top: orbit drawing, solution with (f2,aL) in the left, solution with (Q',^) in the right. The 
dotted curves are the primary and secondary orbits. The stars (*) are at the binary components positions at the time of observations. 
The crosses are the astrometric data. The open squares are the best-fitting model for the solution (0,q:l); the open diamonds for the 
solution (Q',a' L ). Bottom: line visibility (left) and line phase (right) versus spatial frequency radius R in the (u,v) plane. 



Table 5. Retrograde binary scenario fit parameters. 







n = (234 ± 8 


)° 




Q> = (44 ±8)° 






tp,i 


= (1.7 ±0.6) MJD 




= (2.3 ±0.7) MJD 








Xl = 3.2 






Xl = 3.4 






Night 


4>p (°) 


Br7l (A) 


Br72 (A) 


*p (°) 


Br7l (A) 


Br72 (A) 


Br 7 T (A) 


MR1 


255 


-1.5 ±0.4 


-3.2 ±0.9 


244 


-4.0 ±0.7 


-0.8 ±0.1 


-4.8 


MR2 


308 


-0.6 ±0.3 


-3.3 ±2.0 


298 


-3.9 ±0.6 


-0.0 ±0.0 


-3.9 


MR3 


346 


-4.6 ± 1.5 


-2.8 ±0.9 


335 


-4.5 ± 1.3 


-2.9 ±0.8 


-7.4 


HRf 


330 


-5.1 ±0.3 


-4.0 ±0.2 


319 


-5.8 ±0.3 


-3.2 ±0.2 


-9.1 



The table presents, for the two models labeled by Q and the time of periastron passage (tp,i, t' p ., in Modified Julian Date — MJD), 
xl, the position in the orbit (<pp, periastron is 0°) and each component (Br7l, Br72) and total (Br7T) line equivalent width. 

lArtvmowicz fc Lubowl Il996t ) presents a tidally truncated 
disc for a system with a mass ratio 1:1.27 and e = 0.5. 
HD 104237 has both a larger mass ratio (1:1.57) and ec- 
centricity (e = 0.64). Its truncation radius is expected to 
be located sli ghtly further out than in the simulation, at ~ 
0.73 AU using IPichardo et~aTI (|200Sl ) scalings for HD 104237 



parameters. The squared visibility continuum models pre- 
sented in Section [3] have radii smaller than the hydrody- 
namic truncation radius, the same app lying to previous fits 
l|Tatulli et al.ll2007l ; iKraus et al.1 120081 ). They hint at addi- 
tional mate rial in the tidally t r uncated region, s uch as found 
in DQ Tau |Boden et al.ll2009l : ICarr et al)l20oil ). This is fur- 
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Figure 11. Pole-on system diagram in apastron, with horizontal 
arbitrary position angle. The orbit of the primary is in red and of 
the secondary in black. The remaining curves trace, in increasing 
radius order: a) 0.5 AU radius circumbinary disc (magenta); b) 
computed sublimation ring (blue); c) computed sublimation el- 
lipse due to the primary orb i t (red ) . The background is an image 
from [Artvmowicz fc Lubowl jl996h at apocentre. The image was 
scaled-up such that the stars are positioned over the orbits. 



ther supported by the unresolved contribution in the con- 
tinuum models being roughly 50 per cent, l arger than th e 
stellar 30 per cent derived in the SED fits bv lTatulli et all 
In steady state the dust is s ublimated inside a ring 
of radius r BU b ~ 0.6 AU (e.g. iDullemond et all l200ll ; 
llsella fc Nattall2005l ). It is depicted in Fig. [TT] and slightly 
smaller than the tidal t runcation radius. For the tidally dis- 
rupted disc of AK Sco, lAlencar et al.l (|2003l ) also find that 
the sublimation radius is smaller than the truncation radius, 
i.e. the circumbinary disc rim has a temperature below the 
sublimation temperature. Due to the orbital movement of 
the primary, the effective sublimation radius could be larger 
and the dust sublimation "rim" elliptical (cf. Fig. Hip , with 
an eccentricity of 0.36. The previous reasoning is only valid 
if the local cooling time (t C ooi) in the rim i s large r than the 
orbital period (P) . Following iNagel et all (|2010l ) the cool- 
ing time is estimated t coa i ~ 16 d. As the primary orbits, 
the circumbinary disc rim region nearer the star will have 
a higher temperature than the regions further out. This re- 
gion will appear as a hot spot/arc that rotates around the 
disc inner rim, following the primary orbital motion. This 
hot spot/arc could also account for an excess of unresolved 
flux observed in the data. 



6.2 Accretion streamers 



iBeck et al l l|2012l ) detect enhanced H2 emission inside the 
tidally disrupted disc of GG Tau. The peak of H2 emis- 
sion is located at the inner edge of the circumbinary disc 
and coincident w it h a dust accretio n stre amer found by 
iPietu et all l|201ll) . iFolha fc Emerson! (|200ll ) report drastic 
changes in the Br7 emission of GG T au. Streaming activit y 
was also detected in younger systems (|Mavama et aljfeoiot ) . 
High-resolution (magneto-)hydrodynamical models of accre- 



tion streams (e. g. Gunther fc K iev 2002]; Ide Val-Borro et"al] 
l201ll ; IShi et al.ll2012l) find that the streams spiral from the 
circumbinary disc via the Lagrange points into the binary 
components. The mass accretion is variable but the streams 
are stable through the orbits. Our observations find that 
~ 20 per cent of the total continuum is unresolved and not 
accounted by the stellar photospheres. The Br7 data shows 
that > 90 per cent of the emission is compact and there- 
fore only a small fraction < 10 per cent can originate from a 
similar region where H2 is observed in GG Tau or from large- 
scale streams. However, streaming material could contribute 
to the excess continuum and Br7 emission in a shock fur- 
ther down the accretion stream, nearer the stars. The pole- 
on geometry would imply a very small radial velocity for 
the streaming gas and such a contribution (in contrast with 
the jet case) would p ass unnoticed in the total Br7 profile. 
lHanawa et all (|2010h find compression shocks in the region 
of the interaction of each circumstellar disc and transient 
hot spots in the primary circumstellar disc. Such regions 
are also a good locus for compact Br7 emission although, it 
is not obvious if they are directly applicable to HD 104237 
whose stars are expected to be devoid of discs as will be 
argued next. 



6.3 Interacting magnetospheres 

Tidal effects are expect ed to truncate the circ umstellar discs 
from outside. Applving lPichardo et alj l|2005l ) expressions to 
HD 104237, upper limits of 4 R© for the primary disc and 
3 R,0 for the secondary disc are obtained. On the other hand 
the stellar magnetic fields would truncate the circumstellar 
disks from the inside. The search fo r line- of-sig ht magnetic 
fields in HD 104237 bv lDonati et all (Il997l ) and I Wade et al.1 
l|2007l) found marginal detections at a level of 50 G. How- 
ever, a 150 G dipolar field in tensity^ is required t o truncate 
discs of such small radii (e.g. Bess olaz et al.ll2008l ). This dif- 
ference could be reconciled by the system pole-on geom- 
etry, the dipolar magnetospheric loops would be reaching 
the photosphere almost perpendicular to the line-of-sight, 
the detected projected field being small with respect to the 
total dipolar field. Alternatively, there exists a strong field 
but i t is multipolar or not well ordered and of difficult detec- 
tion. iGregorv et all (|2012l ) finds this is indeed the case for 
more massive pre-main-sequence stars with well developed 
radiative cores, such as HD 104237. At the small radii of the 
putative circumstellar disks, multipolar fields can contribute 
significantly to the trunc ation and relax the requir ed dipolar 
field to even lower values (| Adams fc Gregorv||2012h . Circum- 
stellar discs are most probably untenable in each component. 

At periast r on th e distance between the stars is 17 Rq. 
I Get man et al.l (|2008l ) find that the fast rotating disc-less 
young stars have typical flares with sizes similar to the co- 
rotation radius. The size of the magnetospheres can be taken 
as the co-rotation radius. For the primary it is 12 Kq. The 
secondary could be a fast rotator (with a similar co-rotation 
radius as the primary) or a CTTS-like slow rotator with an 
even larger co-rotation radius. The magnetospheres of both 
stars are therefore expected to interact during periastron 



3 Taking iGarcia-Lopez et al.l accretion rate and other relevant 
data in Appendix [Al 
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passage and flaring activity as observed in U ZTau, DQTau, 
AKSco or V4046Sgr would take place (e.g. iGetman et al.l 
120111 ; iKospal et allhoill ; ISalter et all |2010| ). These flares 
would liberate energy and increase the effective size of the 
system magnetosphere. If the Br7 originates in this disor- 
dered magnetosphere, its flux would increase at periastron 
without important changes in line profile, as shown in Sec- 
tion 15.21 and observed in Fig. [4] This B r7 line prof i le be - 
haviour is similar to the one referred to bv lBasri et al.l (|l997l ) 
for the Balmer lines of DQ Tau. 



Table Al. Fit of radial velocity data of lBohm etafl l l2004h . 



Physical quantity 



Value 



t P (HJD - 2 450 000 d) 
K"l (km s" 1 ) 

e 

Wl/2 (°) 
P(d) 

7 (km s 1 ) 



1647.539 ± 0.003 
17.8 ±0.2 
0.643 ± 0.006 
218.2 ±0.9 
19.856 ± 0.002 
14.1 ±0.1 
1.8 



6.4 The Lya jet 

iGradv et al.1 (|2004l ) detect a Lya jet with a speed of ~ 
340 km s _1 in HD 104237. In jet models where the launch 
takes place at the disc surface (cf. iKonigl fc Salmeronll201ll . 
for a recent review) matter is strongly coupled to the mag- 
netic field. For a given field line anchored at a cylindrical 
radius (r^) at the disc surface, most of the acceleration takes 
place up to the Alfven radius (t-a) and is measured by the 
ratio {ta/t-^). If HD 104237 Lya jet is launched from the in- 
ner rim of the circumbinary disc, the ratio can be estimated 
TA/fk — 5. This ratio is also related to the mass loss to mass 
accretion ratio and can be used to estimate the wind mass 
loss M w ~ 2 x 10" 9 M yr" 1 . The magn etic field in the 
launch ing region can be estimated following lAnderson et al.l 
(2005) as Bp ~ 1 G. Such a high field cannot originate in a 
dipole scaling, but is typical of the fields required by mod- 
els for the launching regions of jets in CTTSs. A centrifugal 
driven jet originating in the circumbinary disc is a plausible 
explanation of the Lya jet. 

A disc wind was introduced bv lWeigelt et al.l (|201 lh to 
explain the Br7 of the Herbig Be star MWC 297. But for this 
star the emission is resolved and not compact as we reported 
for HD 104237 in Section [4] Such a model is not directly ap- 
plicable to HD 104237. The absence of significant extended 
Br7 emission could be related to the tidal disruption of the 
circumstellar discs. 

Models that rely on stellar magnetosphere- disc cou- 
pling, such as the X-wind (e. g. ICai et alj l2008l ) or disc- 
magnetosphere simulations (e.g. ILii et all 20121) . are not di- 
rectly applicable to HD 104237 because spiraling material 
crossing the dynamically truncated gap is not expected to 
form stable circumstellar discs (as argued in Section |6.3[) . 
HD 104237 poses a challenge to current disc-magnetosphere 
interaction mo dels of jet s . It w ould be of interest to extend 
the analysis of lShi et al.l (|2012l ) to jet formation. 

The L ya jet cou l d be driven from the star. IUE obser- 
vations by IHu et"ai1 (| 199 ll ) find the Mgll h and k lines to 
have a strong P Cyg profile, with the absorption extending 
up to ~ 520 km s - . This is striking evidence of a powerful 
stellar wind from the polar regions of the primary. Interest 
in stellar winds has been revis ed after the suggestion that 
they tap on accretion energy (|Matt fc Pudritzl l2005l ). De- 
tailed models for winds have been developed and applied 
to slow rotating CTTSs (jCranmedbOplSl . 120081 ) finding mass 
losses lower than observed. However, the HD 104237 primary 
has a higher rotation rate and the expected mass loss would 
increase. On the other hand, the tidally perturbed accretion 
and the magnetospheres interaction at periastron could en- 
hance the amount of energy dissipated in Alfven waves and 
power the wind. An attractive scenario for HD 104237 is one 



consisting of a stellar wind colli mated by an outer c ircumbi- 
nary disc wind, as proposed bv ISautv et al.l (|201ll ) for the 
intermediate mass star RYTau. 



7 CONCLUSIONS 

Spectro-interferometric observations of HD 104237 in the 
Br7 line and adjacent continuum with X/2B — 1.84 mas 
maximum resolution are reported. It is found that the line is 
variable, presenting roughly a factor of two stronger equiva- 
lent width at periastron passage. At least 90% of the the line 
emission is angularly unresolved. The spectro-astrometry of 
the line presents displacements with a position angle com- 
patible with the jet. However, it is argued that emission 
does not originate in the jet but instead in the stars. It is 
found that modeling the spectro-astrometry and the angular 
size of the Br7 with a binary provides good agreement with 
observations. A mixed origin in a binary plus jet was not ad- 
dressed due to the high number of parameters with regard to 
the available data. HD 104237 corroborates that tight bina- 
ries of the Herbig Ae class present similar phenomena as well 
as similar challenges to the disc-magnetosphere interaction 
as their lower luminosity CTTSs counterparts. 
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APPENDIX A: STELLAR PARAMETERS 

The assumed parameters for the HD 104237 system used 
throug hout the paper are presented. A fit to the lBohm et al.l 
(2004) radial velocity data was done to derive error bars on 
the system orbital parameters, the results are presented in 
Table |A"T1 To derive the orbital semi-major axis the masses 
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Table A2. Physical parameters of the system. 



Physical quantity 


Value 


Reference 


mi (Mg) 


2.2 ±0.2 


Bohm et al. 2004 


m 2 (M ) 


1.4 ±0.3 


assumed 


a (AU) 


0.22 ±0.06 


derived 


a (mas) 


1.9 ±0.6 


idem 


in 


1? +12 
i '-9 


idem 


T eff ,i (K) 


8500 ± 150 


Fumel & Bohm 2012 


(L©) 


35t» 


van den Anckcr ct al. 1998 


(R©) 


2.7 ±0.2 


derived 


frot.isini (km s — x ) 


12 ±2 


Donati et al. 1997 


Prot.l (d) 


3.4 ±2.0 


derived 


^C,l/^C,T 


0.2 


Tatulli et al. 2007 




0.1 


ibid. 



Fq i / Fq t is the K-band primary stellar flux to total flux ra- 
tio. Fq 2/ Fc,T IS the secondary to total flux ratio, the remaining 
emission arises in the circumbinary disk. 



must be known. The primary mass is taken as 2 .2 ±0.2 Mp>. 
The mean value is the o ne determined by iBohm et al.l 
and iFumel fc Bohml d2012l). by placing the primary in the 
HR diagram and using IPalla fc Stahierf <|200lh tracks. How- 
ever, the error is increased here from 5 to 10 per cent, 
more in line with discussions of un certainties of pre-main- 
sequence mass determinations (e.g. | Blondel fc Tiin A Diid 
I2006J ; IGennaro et"alll2012l ; iMathieu et al.ll2007l ; ISiessll200 ll ). 
This is still small er than the err or assumed in the more 
recent analysis of IFumel fc Bo hm. which probably overes- 
timates it because of the use of mass-lu minosity rela tions. 
The secondary mass was determined bv lBohm et al.[ to be 
1.7 Mp) , a value quite hig h for a K3 pre-main-sequence star 
(e.g. IMathieu et al.ll2007f ). Furthermore, such a high mass 
implies a system inclina tion angle of 15°. Using the Lya 
jet detection information, iGradv et al.l determine an inclina- 
tion angle of 18 jz^j °. Using pseudo-periodic variability of the 
Ha line and the pr ojected rotational velocity of the system, 
IBohm et al.l |200fj ) derive an inclination of 23^g , which im- 
plies a secon dary mass of ~ 1.2 Mg. On the other hand if 
we place the B5hm et all l|2004h secondary HR position in 
IGennaro et al.l tracks we obtain a mass of < 1.4 M© , high- 
lighting theoretical uncertainties. In this study, a secondary 
mass of 1.4 ± 0.3 Mq is assumed. With the system mass 
fixed, the radial velocity data allows the determination of 
the remaining orbital parameters in Table lA"2l 

For simplicity (and given the available data), full 
coplanarity is postulated, i.e. the stellar, orbital and disc 
inclination angle are assumed to be approximately the 
same. Spectr o -astro metr ic observations of young bi naries by 
iBaines et ail |2006j ) and IWheelwright et alj |201ll ) support 
that the circumprimary disc is coplanar with the binary 
orbit. Observations of more evolved systems find orbital 
and stellar axis align ment for small separations (HaleJ[l994|; 
iHowe fc~C larkc 2009). It is a plausible assumption for a sys- 
tem that probably formed via fragmentation of a common 



APPENDIX B: LEAST SQUARES DISPERSED 
VISIBILITY DECONVOLUTION 

The extraction of the line-to-continuum complex visibility 
from the total to continuum c o mplex visibility was addressed 
previously by IWeigelt et ail (2007), but the effect of the 
spectrograph line spread function was ignored. This is a 
good approximation when the line is well resolved by the 
spectrograph and the SNR is high, which is not the case 
for our observations. In this section we present a procedure 
to 'deconvolve' the complex visibility from the effects of the 
spectrograph and recover the line-to-continuum complex vis- 
ibility. 



Bl Spectral convolution in AMBER data 

The AMBER data consists of an interferometric channel 
and three photometric channels, the four channels being 
spectrally dispersed through a spectrograph. Carrying waves 
are fitted to the interferometric channel, extracting the real 
SH(A) and imaginary 3(A) parts of the coherent flux, at each 
wavelength. Because the data is dispersed, the observables 
are convolved with the spectrograph line spread function 
profile V(X). The interferometric observables are therefore 



{ixThTjVi]} * 7>) ■ r and 3 = ({^/JifjVij} * V) 



where V ftfjVij is the coherent flux vector (in the complex 
plane) for baseline ij and r, i are the versors (in the com- 
plex plane) along the real and imaginary axes. The detected 
photometric signal Fi is also convolved by V 

Fi = fi * V. 

The above quantities are combined in the pipeline to yield 
the squared visibilitjQ 



1 'i. i 



,)*V\ 



The above relation shows that the measured square visibility 
is not the convolution of the visibility by the spectrograph 
line spread function. The delivered differential phase observ- 
able is 



A<t>. 



'T/C 



= arctan|^--tanf^ 



which is also not a simple convolution of a real phase by the 
spectrograph line spread function. To invert the previous 
relations some model for the data has to be assumed. 



B2 Spectral deconvolution method 

The spectra to continuum ratio in the photometric channel 
of telescope i near Br7 is described by 



Fi 



i/C 



l + F L/c P L , 



where Jl/c = Fl/Fc is the ratio of line amplitude to con- 
tinuum value and Pl is the convolution of the intrinsic Br7 



4 Details connected to debi asing the data are ignored as they are 
achromatic in the pipeline jTatulli et al.|[2~007f) . 
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profile pl by the spectrograph line spread function. The co- 
herent flux vector is also assumed to be the combination of 
line and continuum contributions 

\fhflvu = FcVc + F LPL V L . 

The line visibility vector (in the complex plane) Vl can 
change across the profile. However, after verification with the 
HR data, it is assumed to be constant for a given baseline. 
With the above definitions it can be shown that 



Vi 



T/C 



1 + 2F L/C P L V L/C cos(A0 L/r 



) + Fl /c PlV? 



(Bl) 

L/C 



where Vl/c = Vl/Vc, V t/c = Vr/Vc and A0 L /c = 4>^~<t>c- 
The squared visibility normalized by the continuum V^/ C 
is independent of the interferometer transfer function and 
SNR biases and obtained in Section 2. However, we don't 
have direct access to the line-to-cont inuum visibility Vl/c 
and phase A0 L /c jWeigelt et al.ll2007l '). With regards to the 
differential phase 



tan(A<^) T /( 



Vl/cJl/c-Pl sin(A<fc/c) 
1 + Vl/c^l/cPl cos(A0 L /c) ' 



(B2) 



In Eq. IBll and Eg. IB2l the product Fl/c-Pl for a given base- 
line ij can be obtained from the photometric channels data 



Fl/cPl 



(B3) 



For a given baseline, Eqs. IBll to IB3I provide a complete de- 
scription of the data with parameters Vl/c and A0 L /Ci tak- 
ing into account the spectrograph line dispersion function. 



B3 Least squares fitting of the data 



In practice, a Levenberg-Marquardt algorithm is used to fit 
the above model for the V^/ C and tan(A^ T /c) data, taking 
as x = Fl/c-Pl given by Eg. IB3I The x' 2 to be minimized is 



X 



X {V T /c} + X {tan(A0 T/c )} 



The errors in the weighting are constant and estimated 
from the pseudo-continuum rms in the (2.13 — 2.15) ^m 
range. Initial estimates for the parameters are Vl/c — 1 
and A0l/c = 0°, but the results are not sensitive to large 
variations of initial estimates. In the fit only points within 
|xi | ^ 400 km s _1 are used, as including large amounts of 
pseudo-continuum adjacent to the line would bias the re- 
sults. Final xl = 1-3 for the MR data and 1.4 for the HR 
data were obtained. 

The fitted parameters are presented in Fig. [6l the actual 
visibility and phase models are presented in Fig.0and Fig. [3] 
Overall the model is a very good description of the data, 
however there are differential phases (MR1 G1-D0 06:42 and 
MR1 G1-D0 07:18, in Fig. [2) which are slightly red-shifted 
with respect to the line profile and squared visibility. This 
effect is smaller than the error bars in the extracted quanti- 
ties. The data set MR3 H0-G0 07:51 is apparently not well 
fit by the model, it is the worst data in the set with a phase 
rms of 7° . 




Vuc (W) 




A*uc (deg) (W) 

Figure Bl. Top: least squares Vl/c versus differential value 
Vl/c (W). Bottom: least squares A<j!>l/c versus differential value 
A^l/c (W). 



B4 Comparison with the differential method 

The differential method o f lWeieelt et all (|2007T ) was also ap- 
plied to our data. The results of both methods are compared 
in Fig. IBll The differential method presents a result for each 
spectral pixel. In the comparison the average and standard 
deviation of the best points, with \v\ ^ 150 km s _1 , are 
taken. It is found that both methods present the same aver- 
age value within errors. However, the precision of the least 
squares deconvolution is typically a factor of two better than 
that of the differential method. 
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